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Abstract — This paper considers the analysis and optimization 
of a frequency-hopping ad hoc network with a finite number of 
mobiles and finite spatial extent. The mobiles communicate using 
coded continuous-phase frequency-shift keying (CPFSK) modu- 
lation. The performance of the system is a function of the number 
of hopping channels, the rate of the error-correction code, and 
the modulation index used by the CPFSK modulation. For a 
given channel model and density of mobiles, these parameters are 
jointly optimized by maximizing the (modulation-constrained) 
transmission capacity, which is a measure of the spatial spectral 
efficiency of the system. The transmission capacity of the finite 
network is found by using a recent expression for the spatially 
averaged outage probability in the presence of Nakagami fading, 
which is found in closed form in the absence of shadowing and 
can be solved using numerical integration in the presence of 
shadowing. 

I. Introduction 

An ad hoc network or peer-to-peer network comprises 
autonomous systems that communicate without a centralized 
control or assistance. Ad hoc networks include mobile com- 
munication networks that possess no supporting infrastructure, 
fixed or mobile; usually each mobile has identical signal 
processing capability. The preferred channel access for ad 
hoc networks is direct-sequence or frequency-hopping (FH) 
spread spectrum. This paper focuses specifically on frequency- 
hopping spread spectrum ad hoc networks. Such networks 
are characterized by independent, identical, FH radios that 
share the same carriers and frequency channels, and are nearly 
stationary in location over a single hop duration. 

The use of spread spectrum in ad hoc networks provides 
the many advantages of code-division multiple-access com- 
munications |fl~). Two major advantages of FH are that it can 
be implemented over a much larger frequency band than is 
possible with direct-sequence spreading, and that the band 
can be divided into noncontiguous segments. Another major 
advantage is that FH provides resistance to multiple-access 
interference while not requiring power control to prevent the 
near-far problem. Since direct-sequence and narrowband sys- 
tems cannot escape the near-far problem by hopping, accurate 
power control is crucial in limiting the impact of the near- 
far problem in cellular networks [2 |. However, power control 
is not viable for ad hoc networks because of the lack of a 
centralized architecture. Therefore, direct-sequence networks 
require a guard zone imposed by some type of carrier-sense 
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multiple access with collision detection (3j. In contrast, FH 
networks can use an ALOHA-type multiple-access protocol. 

For FH systems, continuous-phase frequency-shift keying 
(CPFSK) is the preferred modulation. Frequency hopping 
with CPFSK offers a constant-envelope signal, a compact 
signal spectrum, and robustness against both partial-band 
and multiple-access interference [2|. CPFSK modulation is 
characterized by its modulation order, which is the number 
of possible tones, and by its modulation index h, which is the 
normalized tone spacing (assumed here to be constant). For 
a fixed modulation order, the selection of h involves a trade- 
off between bandwidth and performance J4). The resistance 
to frequency-hopping multiple-access interference generally 
increases with decreasing h while the error-rate generally 
decreases with increasing h. When the ad hoc network uses 
channel-coded CPFSK modulation, the performance is a func- 
tion of not only the modulation index, but also the rate of 
the error-control code and the number of frequency-hopping 
channels. As each of these parameters influences bandwidth 
and energy efficiency, there is a tradeoff among them, and 
the main goal of this paper is to gain some insight into this 
tradeoff. 

An appropriate objective for the optimization of the net- 
work's operational parameters is to maximize the transmission 
capacity 0, (6), which is a measure of the spatial spectral 
efficiency and is constrained in this paper to explicitly take 
into account the CPFSK modulation format. In contrast with 
the recent literature (e.g., [7], [8|), which assumes an infinite 
number of mobiles spread over an infinite area, in this paper 
the spatial extent of the network and number of mobiles are 
finite. Each mobile has a uniform location distribution with 
an allowance for the mobile's duty factor and shadowing. The 
analysis uses a recent expression [9| for the exact outage 
probability in the presence of Nakagami fading conditioned 
on the network topology, which is determined by both the 
mobile locations and the shadowing factors. The modulation- 
constrained transmission capacity is found by averaging the 
outage probability over the network topoloty, which is possible 
in closed form in certain cases (i.e., no shadowing, and 
interfering mobiles uniformly distributed within an annulus 
centered at the reference receiver) and through numerical 
integration in other cases (e.g., networks with shadowing). By 
maximizing the modulation-constrained transmission capacity, 
the optimal set of operating parameters (code rate, modulation 
index, and number of frequency channels) is found for a given 
channel model and density of mobiles. 



The remainder of the paper is organized as follows. Section 
[iljpresents a system model, which culminates in an expression 
for the instantaneous SINR at the reference receiver. Section 
|ni1 reviews a new expression [9| for the conditional outage 
probability; i.e., the probability that the SINR is below a 
threshold, given a particular set of mobile locations and 



shadowing factors. Section IV discusses how to remove the 



conditioning by averaging over the spatial and shadowing 
distributions. Section |V] discusses the modulation-constrained 
transmission capacity, which is the area spectral efficiency 
of the network and the metric used for the optimization. 



Section VI describes the procedure used for determining the 
combination of parameters that maximize the transmission 
capacity and gives optimization results for several channel and 
network models. 

II. Network Model 

The network comprises M + 2 mobiles that include a 
receiver, a reference (source) transmitter X a , and M in- 
terfering transmitters Xi, Xj^ . The coordinate system is 
selected such that the receiver is at the origin. The variable Xi 
represents both the i th transmitter and its location, and \\Xi\\ 
is the distance from Xi to the receiver. While the interferers 
can be located in any arbitrary region, we assume they are 
located in an annular region with inner radius r ex and outer 
radius r net . A nonzero r ex may be used to model the effects 
of interference-avoidance protocols 0. 

Xi transmits a signal whose average received power in the 
absence of fading and shadowing is Pi at a reference distance 
do- At the receiver, X^'s power is 



(1) 



where gi is the power gain due to fading, £j is a shadowing 
coefficient, and is a path-loss function. Each gi = af, 

where aj is Nakagami with parameter m,-, and E[g$] = 1. 
In Rayleigh fading, m, = 1 and gi is exponential. In the 
presence of log-normal shadowing, the {^} are i.i.d. zero- 
mean Gaussian with standard deviation a s dB. In the absence 
of shadowing, £j = 0. For d > do, the path-loss function is 
expressed as the attenuation power law 



f(d) 



(2) 



where a > 2 is the attenuation power-law exponent, and do 
is sufficiently large that the signals are in the far field. 

Channel access is through a synchronous frequency-hopping 
protocol. The hopping is slow, with multiple symbols per hop, 
which is a more suitable strategy for ad hoc networks than fast 
hopping \2\. An overall frequency band of B Hz is divided 
into L frequency channels, each of bandwidth B /L Hz. The 
transmitters independently select their transmit frequencies 
with equal probability. Let pi denote the probability that 
interferer Xi uses the same frequency as the source. Let 
Di < 1 be the duty factor of the interferer. It follows that 
Pi — Di/L and that using a duty factor less than unity 
is equivalent to hopping over more than L frequencies 0. 



Assuming that Z), = D for all interferers, L' = L/D denotes 
the equivalent number of frequency channels. It is assumed 
that the {gi} remain fixed for the duration of a hop, but 
vary independently from hop to hop (block fading). While 
the {gi} are independent, they are not necessarily identically 
distributed. The channel from each transmitting mobile to the 
reference receiver can have a distinct Nakagami parameter m^. 
The instantaneous SINR at the receiver is 
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Po 
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where Af is the noise power and Ii is a variable that indicates 
the presence and type of interference (i.e. co-channel interfer- 
ence or adjacent-channel interference). When adjacent-channel 
interference is neglected, Ii — 1 when Xi selects the same 
frequency as X , and Zj = otherwise. It follows that Ii is 
Bernoulli with probability P[Ii = 1] = pi. 
Substituting ([]} and Q into the SINR is 

5o^o 



A I 



(4) 
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where T = r^Po/Af is the signal-to-noise ratio (SNR) when 
the transmitter is at unit distance and fading and shadowing 
are absent, fi< = (f , i /P )10 ? ' /10 ||^ir Q! is the normalized 
power of Xi, and fl = 10«°/ 10 | \X Q \ \~ a . 

III. Conditional Outage Probability 

Let /3 denote the minimum SINR required for reliable re- 
ception and ft = {Oq, r^A/} represent the set of normalized 
powers. An outage occurs when the SINR falls below /3. 
Conditioning on ft, the outage probability is 



(o = P[l<P\ft] 



(5) 



Because it is conditioned on ft, the outage probability depends 
on the locations of the mobiles and the shadowing factors, 
which have dynamics over timescales that are much slower 
than the fading. By defining a variable 

M 
i=l 



Z = 



kg&i 



(6) 



the conditional outage probability may be expressed as 

e = Plz^r-^ft] = Fzir-^ft) 



(7) 



which is the cumulative distribution function (cdf) of Z con- 
ditioned on ft and evaluated at T -1 . 

Under the restriction that the Nakagami parameter mo of 
the source's channel is integer-valued, the complementary cdf 
of Z conditioned on ft is shown in |9| to be 
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where /3 = m (3/il , 
% ■■ 

H k (9) = 
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the summation in ( 10 1 is over all sets of positive indices that 
sum to k, and 
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A. Examples 

In the following examples, the source transmitter was placed 
at unit distance from the receiver; i.e., ||Xo|| = 1, and fifty 
interferers were independently placed according to a uniform 
distribution in an annular region with inner dimension r ex = 
0.25 and outer dimension r net = 4. The resulting network 
is shown in the inset of Fig. [T] The $1 was determined by 
assuming a path-loss exponent a = 3, a common transmit 
power Pi = Pq, and the absence of shadowing. The equivalent 
number of frequency channels was set to L' = 200, and the 
SINR threshold set to /3 = 3.7 dB. 

Example #1. Suppose that all signals undergo Rayleigh 
fading. Then mi 



F z (z|0) 



1 for all i, and d8H 1 1 1 specialize to 

JUT ' — * ' ' 

-A,«TT 1 + A)(1 ~ P<)n< 

11 i + AA 



(12) 



The outage probability was found by evaluating (12i at 
z = r _1 and is shown in Fig. [I] Also shown is a curve 
generated by simulation, which involved randomly generating 
the exponentially-distributed {gi} and drawing each inter- 
ferer's hopping frequency from a uniform distribution. The 
analytical and simulation results coincide, which is to be 



expected because ( 12 1 is exact. Any discrepancy between the 
curves can be attributed to the finite number of Monte Carlo 
trials (one million trials were executed per SNR point). 

Example #2. Now suppose that the link between the source 
and receiver undergoes Nakagami fading with parameter too = 
4. In this case, (|8|l-(fTT|> must be used to find the conditional 
outage probability. Using the same (3, L', and f2 as Example 
#1, the outage probability was found and plotted in Fig. 
[T] Two cases were considered for the interferer's Nakagami 
coefficient: = 1 and vrn — 4. The mi — 4 case represents 
the situation where the source and interferer are equally visible 
to the receiver, while the m; = 1 case represents a more typical 
situation where the interferers are not in the line-of-sight. As 
with the previous example, the analytical curves are verified 
by simulations involving one million Monte Carlo trials per 
SNR point. 
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Fig. 1. Conditional outage probability e of Examples #1 and #2 as a function 
of SNR r. Analytical curves are solid while dots represent simulated values. 
Top curve: all channels Rayleigh. Middle curve: all channels m = 4. Bottom 
curve: m = 4 from source, and m = 1 from interferers. The network topology 
is shown in the inset. The receiver is represented by the star at the center of 
the radius-4 circle while the 50 interferers are shown as dots. 

IV. Spatially Averaged Outage Probability 

Because it is conditioned on J~2, the outage probability 
e Q presented in the last section depends on the location of 
the interferers and the values of the shadowing factors. The 
conditioning on J~2 can be removed by averaging Fz(z\tt) with 
respect to fi, which results in the average outage probability 

e = E n [e ] = E[F z (T- 1 \Si)]=Fz(T- 1 ). (13) 

Finding Fz {z) = 1 — F Z (z) requires the following integration: 

Fz(z) = [ (n/n>i)Vz(*Hdw (14) 



\i=0 



where famuli) is the pdf of fij and the {fli} are assumed to 
be independent. 

In (9), it is shown that F z (z) in the absence of shadowing 
is given by ( 15 1 at the bottom of the page, where q = Pq/Pi, 



J{y) = 2F1 
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and 2-F1 is the Gauss hypergeometric function, which has the 
integral representation 



2 Fi([a,b];c;x) = 
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While the hypergeometric function is itself an integral, it is 
widely known and is implemented as a single function call in 
most mathematical programming languages, including Matlab. 

When shadowing is present, Q$ = 10^*' 10 ||Xi||~ 
1, which has pdf 

[C (c^C t )-C(c^r«)] 
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for < bj < oo, and zero elsewhere, where 

.2/ 



CO) = erf 



aj ln^(10) - 50a In (z) 
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Because ||X || is deterministic, D, Q = 10« o/10 | |X | |~ Q is a 
log-normal variable with pdf 
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ln(10)w 

for < lj < oo, and zero elsewhere. 

By substituting p0| ) and ( fT8j ) into ( fl~4| ) and using the 
definition of (3 , the complementary outage probability is 
found to be (f22j), given at the bottom of the page, where 
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In order to compute the outage probability, the integral 



inside the product in (22i can be evaluated numerically by 



Simpson's method, which provides a good tradeoff between 
accuracy and speed, while the second integral can be evaluated 
through Monte Carlo simulation. 

Example #3. The average outage probability is shown as 
a function of M for two L' in Fig. [2] The M interferers are 
uniformly located within the annulus bounded by r ex = 0.25 
and r net = 4. The results assume a path-loss exponent a = 3, 
a common transmit power P t = P , and that the source is at 
unit distance. A mixed-fading model is assumed, i.e. tuq = 4 
and m l = 1 for i > 1, the SNR is T = 10 dB, and the SINR 
threshold is j3 = 3.7 dB. For each value of L', four curves are 
shown, one without shadowing and three with shadowing. To 
model log-normal shadowing, each £j drawn from a zero-mean 
Gaussian distribution with standard deviation a s that varies 
among the set a s = {2,4,8} dB. From Fig. [2] it is observed 
that the outage probability degrades with increasing M and 
decreasing L', and that shadowing increases the average e 
increasing <j s for the given parameters. 
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Fig. 2. Average outage probability as a function of the number of interferers 
M for two values L' . For each L', curves are shown for the case of no 
shadowing, and for shadowing with three values of <j 3 . 

V. Transmission Capacity 

Often, networks are constrained to ensure that the outage 
probability e does not exceed a maximum outage probability ( 
£ [0, 1]; i.e., e(A) < £ with the dependence of e on the density 
of transmissions A made explicit. Under such a constraint, 
the maximum density of transmissions is of interest, which 
is quantified by the transmission capacity (TC) Q, O. With 
outage constraint (, the TC is 



(0 



(25) 



where e _1 (C) is the density of the underlying process whose 
spatially averaged outage probability satisfies the constraint 
e(A) < C Wim equality, and (1— C) ensures that only successful 
transmissions are counted. The TC is a measure of the spatial 
intensity of transmissions and has the units of number of 
(successful) transmissions per unit area. 

As originally defined in |5|, the transmission capacity is 
a function of the SINR threshold j3 and is found without 
making any assumptions about the existence of any particular 
type of modulation or channel coding. In practice, the SINR 
threshold is a function of the modulation and coding that is 
used. Let C(j) denote the maximum achievable rate that can 
be supported by the chosen modulation at an instantaneous 
SINR of 7 assuming equally likely input symbols; i.e., it is 
the modulation-constrained capacity or symmetric -information 
rate. If a rate R code is used, then an outage will occur 
when C(7) < R. Since C(^) is monotonic, it follows that 
j3 is the value for which C((3) — R, and therefore we can 
write f3 = C~ 1 (R). If it is assumed that each codeword is 
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transmitted using one hop, then the outage probability with 
SINR threshold j3 = C^ 1 (R) provides a good prediction of 
the codeword error rate. 

Frequency-hopping systems often use noncoherent CPFSK 
modulation Ipfl, El, The maximum achievable rate of nonco- 
herent CPFSK is given in |4] for various modulation indices h, 
where it is called the symmetric information rate. In particular, 
Fig. 1 of |4| shows the symmetric information rate of binary 
CPFSK as a function of 7 for various h. To emphasize the 
dependence of the capacity on h, we use C(h, 7) in the sequel 
to denote the rate of CPFSK with modulation index h. For 
any value of h, the value of the SINR threshold j3 can be 
found from the corresponding curve by finding the value of 
7 for which C(h, 7) = R. For instance, when R — 1/2 
and h = 1, the required /3 = 3.7 dB. In [10|, it was found 
that in practice, and over a wide range of code rates, turbo- 
coded noncoherent CPFSK is consistently about 1 dB away 
from the corresponding modulation-constrained capacity limit. 
Thus, the j3 required in practice will generally be higher than 
the value obtained from the capacity limit by a small margin. 
For instance, if a 1 dB margin is used, then the SINR threshold 
for noncoherent binary CPFSK with R = 1/2 and h = 1 
should be set to (3 = 4.7 dB. 

When accounting for modulation and coding, the maximum 
data transmission rate is determined by the bandwidth B /L of 
a frequency channel, the spectral efficiency of the modulation, 
and the code rate. Let 77 be the spectral efficiency of the 
modulation, given in symbols per second per Hz, and defined 
by the symbol rate divided by the 99 percent-power bandwidth 
of the modulation^] The spectral efficiency of CPFSK can 
be found by numerically integrating the normalized power- 
spectral densities given in EJ, or since we assume many sym- 
bols per hop, by Equation (3.4-61) of 1111 and then inverting 
the result. To emphasize the dependence of 77 on h, we denote 
the spectral efficiency of CPFSK as r/(h) in the sequel. When 
combined with a rate-i? code, the spectral efficiency of CPFSK 
becomes RrjQi) (information) bits per second per Hz, where 
R is the ratio of information bits to code symbols. The data 
rate supported by the channel is Ri](h)B/L bits per second. 
The average data rate, or throughput, must account for the 
duty factor D and only count correct transmissions. Hence, 
the throughput is 
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(26) 



L L' 
The modulation-constrained transmission capacity is the 
throughput multiplied by the node density, 

\Rr}{h)B(l - e) 



r(A) = XT = 



(27) 



rg/)) is the number of interferers 



where A = M/(ir(rl 
per unit area. In contrast with EH) , this form of transmission 
capacity explicitly takes into account the code rate R, as 
well as the spectral efficiency of the modulation rj(h). It 

'Percent-power bandwidths other than 99 can be used, but will influence 
the amount of adjacent-channel interference. 



furthermore accounts for the hopping bandwidth B/L'. Rather 
than constraining outage probability, r(A) and the outage prob- 
ability vary with the node density A. Since it accounts for the 
actual system bandwidth B, (27 1 assumes units of bps/m 2 . By 



dividing by bandwidth, the normalized modulation-constrained 
transmission capacity 



r \Rr,(h)(l 
r(A) = B = LJ- 
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takes on units of bps/ Hz/m . However, unlike (25 1, r'(A) is 



in terms of information bits rather than channel bits. 

VI. Network Optimization 

For a particular spatial distribution and channel model, there 
will be a set of {V 1 R, h) that maximizes the normalized TC 
t'(A). Finding the optimal set of parameters is challenging 
because the search space is large. However, in [12], we 
found that the optimization surface is convex in the case of 
unshadowed Rayleigh fading, and further experiments have 
confirmed that the surface is still convex in the presence of 
shadowing and Nakagami fading. 

When the search space is convex, the optimal set of 
(L', R, h) can be found as follows: 

1) For each L', R, and h select the endpoints for the search 
interval. Initially, the endpoints are widely separated to 
ensure that the optimal point lies in the search interval. 

2) Compute r'(A) for the endpoints and the midpoint of the 
interval of one of the parameters in the set (L',R,h). 

3) For the same parameter of the set (L',R,h) used in 
step [2] move the midpoint of the interval towards the 
endpoint that gives higher r'(A). Reduce the size of the 
interval by moving the endpoints closer. 

4) Repeat step [2] and [5] recursively for all the parameter in 
the set (L' , R, h), until the following conditions are both 
satisfied: 

a) The midpoints of all three intervals are the same 
as in the previous iteration; 

b) The difference between the endpoints and the mid- 
point of the interval is equal to a given tolerance. 

Optimization results were obtained for M — 50 interferers 
at a fixed SNR of Y = 10 dB. The networks had an inner radius 
of r ex = 0.25 and one of two maximum radii: r net = 2, which 
creates a dense network, or r net = 4, which creates a sparser 
network. Three sets of fading coefficients were considered: 
m» = 1 for all transmitters (universal Rayleigh fading), 
rrii — 4 for all transmitters (universal Nakagami fading), and 
a combination of too = 4 for the source and rrii = 1 for 
the interferers (mixed fading). Unshadowed (a s = dB) and 
shadowed (<j s — 8 dB) environments were examined. 

The results of the optimization are shown in Table [j] For 
each set of parameters, the (L 1 , R, h) that maximize the 
TC are listed, along with the corresponding normalized TC 
r' t . In addition to showing the TC when using the optimal 
parameters, the normalized TC r' sub is shown for a typical 
choice of parameters: (L', R, h) = (200, 1/2, 1). 



TABLE I 

Results of the Optimization for M = 50 interferers. The 

NORMALIZED TC T 1 IS IN UNITS OF BPS/KHZ-m 2 . 
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When R is related to j3 through the modulation-constrained 
AWGN capacity, the optimization assumes Gaussian interfer- 
ence and capacity-achieving codes. However, as mentioned in 
Section [V] actual systems require a threshold that is higher 
than this ideal value by some margin. The column labeled t[ 
gives the normalized TC when a 1 dB margin is used with the 
listed values of (L',R,h). A modest loss in TC is observed 
when this margin is used. 

The results shown in Table [I] highlight the importance of 
parameter optimization. The TC is improved by a factor of 
5-10 by selecting optimal, rather than arbitrary parameters. 
The influence of the fading distribution can be observed. 
Performance is worst under the assumption that all channels 
undergo Rayleigh fading, but this is a pessimistic assumption. 
Performance improves when all channels have a common 
m, = 4, but it is unrealistic to assume that all interferers 
are in the line-of-sight. Performance is further improved when 
the interferers undergo Rayleigh fading while the source 
undergoes better than Rayleigh fading (mo = 4). 

The optimal modulation index is h = 0.59 for all network 
configurations, which coincides with the optimal value found 
in I0] for Rayleigh fading. The optimal L' is larger for denser 
networks and when the signals become more line-of-sight. The 
denser network requires a higher R than the sparser network. 
At first inspection, the last result may seem counter-intuitive. 
One would ordinarily expect a denser network to require the 
error-protection of a lower-rate code. However, the sparser 
network uses a lower value of L', and hence, has a higher 
probability of collision pi = 1/L'. The lower-rate code used 
by the sparse network helps to offset this higher collision 
probability. 

Table [I] shows that performance actually improves in the 
presence of shadowing. This is because with shadowing, the 
outage probability is less sensitive to the choice of code rate 
R and the number of hopping channels L'. While the outage 
probability is worse with shadowing than without for a wide 
range of R and L', when the value of R is sufficiently high 
and/or the value of V sufficiently low, the outage probability 
is improved with shadowing. As can be seen in Table [I] the 
optimal values of R and L' are larger and smaller, respectively, 
with shadowing than without. As indicated by ( |27) , increasing 
R or decreasing L' will increase the transmission capacity for 



a given outage probability and modulation index. 

VII. Conclusion 

Frequency-hopping ad hoc networks play an important role 
in modern communications. When used with coded CPFSK, 
which is the typical choice of modulation, the performance 
depends on the modulation index, the code rate, and the 
number of frequency-hopping channels. These parameters are 
often chosen arbitrarily. The procedure disclosed in this paper 
enables the optimization of the parameters for a fixed network 
with a fixed number of users in the presence of shadowing 
and Nakagami fading. The key innovation facilitating the 
optimization is a new closed-form expression for the outage 
probability in the presence of Nakagami fading. 

The results presented in this paper are just a sample of 
what is possible using this methodology. For instance, other 
types of modulation and reception could be considered, such 
as nonbinary CPFSK or multi-symbol noncoherenet reception 
fl3l . The influence of the spatial model may be studied, as can 
the use of directional antennas. The role of adjacent-channel 
interference can be studied by considering the effect of spectral 
splatter. The influence of receiver position can be studied by 
allowing it to move radially from the center of the network to 
the outer perimeter. 
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